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ABSTRACT 

Recent hydrodynamic studies of core-collapse supernovae imply that the neutrino-heated ejecta 
from a nascent neutron star develops to supersonic outflows. These supersonic winds are influenced 
by the reverse shock from the preceding supernova ejecta, forming the wind termination shock. We 
investigate the effects of the termination shock in neutrino-driven winds and its roll on the r-process. 
Supersonic outflows arc calculated with a semi-analytic neutrino-driven wind model. Subsequent 
termination-shocked, subsonic outflows are obtained by applying the Rankinc-Hugoniot relations. We 
flnd a couple of effects that can be relevant for the r-process. First is the sudden slowdown of the 
temperature decrease by the wind termination. Second is the entropy jump by termination-shock 
heating, up to several 100 A^a^- Nucleosynthesis calculations in the obtained winds are performed to 
examine these effects on the r-process. We find that 1) the slowdown of the temperature decrease plays 
a decisive roll to determine the r-process abundance curves. This is due to the strong dependences of 
the nucleosynthetic path on the temperature during the r-process freezeout phase. Our results suggest 
that only the termination-shocked winds with relatively small shock radii (~ 500 km) are relevant for 
the bulk of the solar r-process abundances {A « 100 — 180). The heaviest part in the solar r-process 
curve {A 180 — 200), however, can be reproduced both in shocked and unshocked winds. These 
results may help to constrain the mass range of supernova progenitors relevant for the r-process. We 
flnd, on the other hand, 2) negligible roles of the entropy jump on the r-process. This is a consequence 
that the sizable entropy increase takes place only at a large shock radius (> 10, 000 km) where the 
r-process has already ceased. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — stars: neu- 
tron — supernovae: general 



1. INTRODUCTION 

Whether the rapid-neutron capture (r-process) can 
take place under the circumstance of core-collapse super- 
novae has been an unresolved problem for a long time (see 
Cowan fc ThielemannI 12004 IWanaio fc Ishimarul l2006l : 
Arnould et al.l I2007L for a recent review). The req- 
uisite physical conditions include neutron-rich environ- 
ment, high entropy, and short expansion timescale to 
obtain the high enough neutron-to-seed abundance ra- 
tio (~ 100) for the production of heavy r-process nu- 
clei. The neutrino-heated supernova ejecta (neutrino- 
driven winds) from a proto-neutron star has been ex- 
pected as a suitable astrophysical site to fulfll such 
physical conditions (in pa rticular high entropy and 
short expan sion time scale. IWooslev fc Hoffma n] |1992l: 
Mever et all [199 2: Wo oslev et al.1 1199 4: Takahas hTet all 
1994 iQia n fc Wooslc/ [l996l: [Card aU fc Fuller HWA 
Otsuki et a l. 2000; Su mivoshi et al.l 12000: Wanaj o et all 
20011: [Thompson et al.ll2001f) . 

Previous studies confirmed, however, that the spher- 
ically expanding winds from a typical proto-neutron 
star (e.g., 1.4 Mq with a 10 km radius) cannot 
attain the requisite physical conditions for the r- 
process. Some mechanisms to cure this problem 
have been proposed, e.g., general relativistic effects 



([Burrows et al.ll2007[). however, no cons ensus has been 
achieved. Recentlv. lArcones et al.l ([20071 ) have suggested 
that the reverse shock propagating from the preceding 
supernova ejecta might have important effects on the 
r-process nucleosynthesis. Two-dimensional hydrody- 
namic studies of core-collapse supernovae have shown 
that the outfiow driven by neutrino heating develops 
to be supersonic, which eventually dece lerate by the re- 
yerse shock f rom the outer layers ( e.g.. iJanka fc Miilleil 
1995, 19 961: [Burrows et all Il995t iBuras et a l. 2006). 

( 2007t l have explored the effects of the re- 
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(iCardafl fc Fullej[l99^lOtsuki et al.l[20 00l: IWanaio et al 
2001'; 'Thompson et al.|l2001'), magnetic field ([Thompson 
2003; Suzuki fc Nagataki 2005), anisotropic neutrino 
radiation ([Wanaiol |2006bj), and acoustic heating 



verse shock on the properties of neutrino-driven winds by 
one-dimensional, long-time hydrodynamic simulations of 
core-collapse supernovae. They have found that, in their 
all models (10 — 25 Mq progenitors), the outflows become 
supersonic and form the "termination shock" when col- 
liding with the slower preceding supernova ejecta. This 
condition continues until the end of their computations 
(10 seconds after core bounce) in their all "standard" 
models (with reasonable parameter choices). Their re- 
sults show that the entropy abruptly increases and the 
temperature continues to d ecrease very s l owly i n the 
termination-shocked winds. lArcones et al.l ([2007[ ) have 
suggested that the strongly time-dependent behaviours 
of the shocked conditions would play a non-negligible role 
to the r-process nucleosynthesis. 

Systematic nucleosynthesis calculations with the 
termination-shocked, supersonic neutrino-driven winds, 
however, have been absent, despite a number of works 
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with t he subsonic ( "bree z e" ) ou tflows fe.g.. lWanaio et all 
1200 It ). iThompson et all ()2001h briefly discussed possible 
effects of the shock-decelerated winds on the r-process 
abundances. They calculated, however, only one spe- 
cific case with an arbitrary chosen shock radii, in which 
the entropy increas e wa s only modera te (~ 10-/VaA:). 
IWanaio et all (|2002f) and lWanaiol (|2007f ) explored the ef- 
fects of the terminated supersonic winds on the r-process 
by introducing a constant "freezeout temperature" Tf. 
They showed the impact of Tf on the r-process abun- 
dance curves, in particular, near the third r-process peak 
(A « 195) and beyond. However, T{ was taken to be an 
arbitrary chosen free parameter. Their approach only 
roughly mimicked the wind-termination effects, and the 
resulting entropy increase was not taken into account. 

In this study, we investigate the effects of the wind 
termination shock on the r-process nucleosynthesis in 
some detail. For this purpose, the thermodynamic 
histories of supersonic outflows are calculated with a 
semi-analytic, spherically syr nmetric, gen eral rclativistic 
neutrino-driven wind model (jWanaio et al . 2001, 200^). 
The termination-shocked outflows are then obtained by 
applying the Rankine-Hugoniot relations at arbitrary 
chosen shock radii (§ 2). In § 3, the behaviors of the 
shock-heated winds are presented along with their de- 
pendencies on the neutron star masses, the neutrino lu- 
minosities, and the shock radii. By applying the ob- 
tained wind trajectories, the r-process calculations are 
performed with the extensive nuclear reaction network 
(§ 4). We discuss effects of the termination-shock decel- 
eration and heating on the r-process nucleoynthesis. Our 
conclusions of this study and discussion are presented in 
§5. 

2. NEUTRINO-DRIVEN WIND MODEL 
2.1. Wind Equations 

After the gravitational collapse of a massive star (> 
8Mq), a great amount of neutrinos are emitted from 
the proto-neutron star. Beyond the neutrino-sphere, 
where the neutrino opacity becomes < 1, heating of 
matter by neutrinos dominates over the cooling process. 
Matter in the vicinity of the neutrino sphere is then 
blown off as "neutrino-driven winds" . The physical pro- 
files of each wind can be obtai n ed b y solving the wind 
equations ( Dunca n et al.l Il986l: iQian & Wooslcv Il996t 
Cardal l & Full er„1997l: |Otsuki et aLf l200Q: ,Wanaio et all 
20011: IThomps on et al.ll2001h . 

Steady state spherically symmetric wind equations in- 
cluding general rclativistic effects are described by the 
following mass, momentum and energy conservations, 
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where r, p, P, e, M, and Mns denote the distance from 
the center, matter density, pressure, specific internal en- 
ergy, mass ejection rate from the surface, and gravita- 
tional mass of the proto-ncutron star, respectively. The 
radial velocity u is related to the proper velocity of the 



matter v as measured by a local, stationary observer by 

+ {u/cf ~ 2GMNs/rc^. The specific heat- 
ing rate q includes both heating and cooling by neutrino 
interactions described in Otsuki et al. (2000), in which 
general rclativistic effects are explicitly taken into ac- 
count. In this study, the neutrino luminosities L^, of all 
flavors are assumed to be equal, and the rms average 
neutrino energies are taken to be 10, 20, and 30 MeV for 
electron, anti-electron, and the other flavors of neutrinos, 
respectively. 

To close the above equations (l)-(3), equations of state 
have to be added. We assume that the wind matter is 
composed of rclativistic electrons, positrons and photons, 
and also non-relativistic free nuclcons. The equations of 
state can be then described as 



P 



180 c3F~ 



11^2 fc4 

60 c^F 



rp4: 



pkT 

f — , 

TUN 

3 kT 

2 TOjv ' 



(4) 
(5) 



where mjv is the nucleon rest mass and T is the matter 
temperature. Effects of arbitrary relativity and degener- 
acy of electrons are ignored, which do not significantly 
modify the thermodynamic histories of winds (e.g., up to 
10% in asymptotic entropy, Wanajo et al. 2001). Solv- 
ing equations (l)-(5) with boundary conditions, we can 
obtain the thermodynamic trajectories 'u(r), T(r), and 
p{r). 

The radius of the proto-neutron star i?NS is taken to 
be 10 km, which is assumed to be equal to the neutrino 
sphere. The boundary conditions at i?NS £^re set to be 
p{Rns) = lO^*' g cm^'^ and T(i?Ns) that equalizes the 
neutrino heating rate with cooling rate. The velocity 
u(i?Ns) (or equivalently Al in eq. [1]) is determined by 
an iterative relaxation method with equations (l)-(5) and 
the above boundary conditions to obtain the transonic 
(i.e., becoming supersonic through the sonic point) wind 
solutions. 

2.2. Termination- Shock Conditions 

In each transonic wind, we assume that a termination 
shock appears at an arbitrarily radius Rs beyond the 
sonic radius Rc, at which the Rankine-Hugoniot shock 
jump conditions are applied. The Rankine-Hugoniot 
equations for mass, momentum, and energy conserva- 
tions are 

Piui = po^o, (6) 
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where the subscripts 1/0 denote quantities in 
shocked/unshocked outflows just above/behind Rs- 
We neglect the rclativistic effects (i.e., u = v) because of 
the distant shock positions {Rs > Rc) from the surface, 
where u < 0.1c (see Figs. 1 and 2). In equations (6)-(8), 
the velocity of the shock position is assumed to be neg- 
ligible compared to the wind velocity lO^cms"^). 
The termination-shocked subsonic solutions are then cal- 
culated by solving equations (l)-(5) with the boundary 
conditions just above Rg obtained from equations (6)- (8) 
(with eqs. [4] and [5]). 
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3. PROPERTIES OF TERMINATION-SHOCKED WINDS 

The solutions of termination-shocked winds u [top 
left), T {top right), and p (bottom left), together with 
the entropy per baryon s {bottom right), as functions of 
r are displayed in Figure 1. The neutron star mass is 
taken to be Mns = 1.4 Mq. The neutrino luminosities 
are assumed to be Li^si = 10 {blue lines) and 1 {green 
lines), where = L^/10^^ erg. The former and lat- 

ter can be regarded as representative of the early and 
late wi nds, respectively (~ 1 — 2 s and > 10 s af ter core 
bounce. IWooslev et al.l ll994':'Arcones et al ."2007^1. In the 
hydrodynamic results by Arcones ct al. (2007), the ter- 
mination shock appears at — 3000 — 7000 km in their 
fiducial model (for the 15 M0 progenitor star), but Rs is 
highly progenitor dependent as well as time dependent, 
ranging from a few 100 km to several 10000 km. We con- 
sider, therefore, various cases with R^ — 500, 1000, 3000, 
10000, and 30000 km, whenever Rg > Rc is satisfied. For 
comparison, the cases with Rs = 00 and Rg = -Rc are also 
displayed in Figure 1 (sonic points are denoted by stars). 
The former corresponds to the unshocked transonic so- 
lution. In the latter case, the shock jump disappears 
and thus the wind is smoothly connected to the subsonic 
breeze solution through the sonic point, with the maxi- 
mum u(i?Ns) (i.e., maximum M). 

As seen in Figure 1, supersonic winds abruptly decel- 
erate at Rs {top left) to become subsonic breeze flows. 
As a consequence, the temperature {top right) and den- 
sity {bottm left) decrease rather slowly after the shock 
jumps, and the entropy increases by the termination- 
shock heating. The shock jump becomes large as the 
shock position moves away from the sonic point. In par- 
ticular, the extreme entropy jump is seen in the distant 
Rs {> 10000 km) wind, up to ~ 400 in the high- 
est case. At a fixed Rs, the early wind (ii/,51 — 10) 
results in a larger entropy jump, which has the larger 
kinetic energy to be converted to thermal energy (Fig. 1; 
top left). The termination-shocked entropies in the early 
winds are appreciably higher, despite their significantly 
lower unshocked entropies than t hose in the late winds. 

The hydrodynamic study by lArcones et al.l ()2007D 
shows variations of M^jg with time as well as depend- 
ing on the progenitor masses, ranging from 1.1 Mq to 
2.0 Mq. In Figure 2, the wind solutions for Mns = 
2.0 M0 {red lines) as an extreme case are compared to 
those for A/ns — 1.4 M© {green lines), in which the neu- 
trino luminosity is fixed to L^^^i — 10 (otherwise the 
same as Fig. 1). 

As can be seen, the unshocked asymptotic entropies 
for Mns = 2.OM0 are about a factor of two higher 
than those for Mns = 1.4 M© (Fig. 2, bottom left) ow- 
ing to general relativistic effects ( Cardall & Fuller 1997; 
Otsuki et al.ll2000HWanaio et al.ll2001l : iThompson et all 
20011) . In addition, we find larger entropy jumps for 
Mns = 2.0 M0 because of the somewhat faster radial 
velocities (Fig. 2, top left). As a consequence, the 
termination-shocked entropies for Mns — 2.0 Mq are ap- 
preciably higher than those for Mns = 1-4 Mq. 

All the above aspects are qualitatively s imilar to those 
found in the hydrodynamic results by lArcones et al.l 
(|2007D . For example, in their fiducial model (with Mns = 
1.2 — 1.3 Mq resulting from the 15 Mq progenitor), the 
unshocked and shocked entropies at 2.0 s after core 



bounce (Li.,51 ~ 10 and Rs ~ 4000 km) are ~ 70 
and ^ 180 NAk, respectively. This is in reasonable agree- 
ment with our result in the wind with Mns = 1.4 Mq, 
= 10, and Rs = 5000 km (not displayed in Figs. 1 
and 2 but shown in Fig. 3), in which the unshocked and 
shocked entropies are 76 NAk and 166 A^a^, respectively. 
Agreements to a similar extent can be found in other 
relevant cases. Our semi-analytic approach is therefore 
quite useful to obtain the thermodynamic histories of 
termination-shocked winds, which are consistent with de- 
tailed hydrodynamic results at qualitative levels. 

We note that the termination-shocked outfiow merges 
with the preceding supernova ejecta and it cannot be re- 
garded as a wind any more. In fact, the hydrodynamic 
calculations show that the shocked outflows further de- 
celerate to be nearl y constant velocities (e.g.. Fig. 5 in 
lArcones et al.l [20071) when accumulated in a dense shell 
between the forward and reverse shocks. Nevertheless, 
we apply the wind solutions for the termination-shocked 
outflows as described in § 2. As a consequence, the out- 
flow velocity rapidly drops (Figs. 1 and 2), resulting in 
the slower decreases of temperature and density with 
time (Fig. 5) t han those in the hy drodynamic calcula- 
tions (Fig. 6 in lArcones et al.|[2007l ). It should be noted, 
however, that the r-process nucleosynthesis ceases only 
within ^ 0.2 — 0.3 s after (or, in some cases, even before) 
the wind termination in the current nu cleosynthesis cal- 
culations (§ 4, see also IWanaioll2007l ). Our treatment 
may be thus reasonable for the current purpose (i.e., for 
nucleosynthesis), since the time variations of tempera- 
ture and density during such a short timescale would be 
negligible in the shock-decelerated outflows. 

Figure 3 shows relations between the entropies Sg and 
the temperatures Tg^g just above the shock radii for var- 
ious Rs cases. The neutron star masses are taken to be 
Mns = 1.4 Mq {crosses) and 2.0 Mq {asterisks). The 
results with various Rs (same as in Figs. 1 and 2) are 
connected with solid and dot-dashed lines for L^^^i = 10 
and 1, respectively. The cases of Rs = -Rc, 1000 km, and 
10000 km are marked by circles, squares, and triangles, 
and the numbers indicated at the circles are the sonic 
radii -Rc. Note that the entropy just behind the shock 
for each case is the same as Sg in the Rs = Rc (sub- 
sonic) wind (marked by circles). As seen in Figure 3, 
the entropy jump occurs only when the temperature de- 
creases below Tg « 2.5 (at which the r-process begins). 
We find a trend that Tg^g is systematically higher in ear- 
lier (i.e., higher L^) winds, while A/ns is less sensitive 
to Tg^s- As an extreme case, we also display the results 
with Mns = 1.4 Mq, -Rns = 10 km, and L^^^i = 70 
{dashed line). In these cases, the entropy jump occurs at 
a smaller Rs and thus at significantly higher temperature, 
but the resulting entropy is only moderate (< 100 NAk) 
for Tg s > 2.5. In addition, the neutrino sphere at earlier 
times would be significantly larger than 10 km and thus 
the entropy would be lower than the current case. 

Moreover, matter in the earlier winds with -L^ si > 10 
would be proton-rich (jBuras et al.l 120061 : lArcones et al.l 
|2007() . in which no r- processing is expected. 

We conclude, therefore, the entropy jump by 
termination-shock heating does not help to enhance the 
neutron-to-seed ratio prior to the r-process (Tg > 2.5). 
The high entropy should be attained much earlier than 
the current situation to infiuence the neutron-to-seed 



4 



Kuroda 



ratio (iQian & WooslevI 119961: 


Suzuki & Nagatakil 120051: 


IWanaiol l2006bl: iBurrows et al. 


I2007D. The termination 



shock can, however, play a decisive role to determin the 
r-process curve as described in § 4. 



4. NUCLEOSYNTHESIS 
4.1. Nuclear Reaction Network 

Adopting the thermodynamic trajectories discussed in 
§ 2 for the physical conditions, the nucleosynthetic yields 
are obtained by solving an extensive nuclear reaction net- 
work code. The network consists of 6300 species between 
the proton and neutron drip lines predict ed by a recent 
fully m icroscopic mass formula (HFB-9, iGorielv et al.l 
[2005a|), all the way from single neutrons and protons 
up to the Z ~ 110 isotopes. All relevant reactions, 
i.e. (n,7), (p,7), (a,7)> (P,"), {a,n), {a,p), and their 
inverse are included. The experimental data whenever 
available and the theoretical predictions for light nu- 
clei {Z < 10) are taken from the REACLIB^ com- 
pilation. The three-body reaction a(an,7)^Be, which 
is of special importance as the bottleneck reaction to 
heavier nuclei, is t aken from the experimental data of 
lUtsunomiva et aL ( 2001). All other reaction rates are 
taken from the Ha user-Feshbach rates of BRUSLIB^ 
( Aikawa et al.ll2005f ) making use of experimental masses 
( Audi et ani2003f) wheneve r availa ble or the HFB-9 mass 
predictions (jGorielv et al.l l2005a[ ) otherwise. The /3- 
decay rates are taken from t he gross theory predictions 
(GT2. lTachibana et al.l[l990l) . obtained with the HFB-9 
Qi3 predictions (T. Tachibana 2005, private communica- 
tion). Neutrino-induced reactions as well as the nuclear 
fission are not considered in this study. 

Each calculation is initiated when the temperature de- 
creases to Tg = 9 (where Tg = T/IO^K). At this high 
temperature, the compositions in the nuclear statistical 
equilibrium are obtained (mostly free nucleons and a par- 
ticles) immediately after the calculation starts. The ini- 
tial compositions are thus given by X„ — 1 — Y,, and 
Xp = Ye, respectively, where Xn and Xp are the mass 
fractions of neutrons and protons, and Kg is the ini- 
tial electron fraction (number of proton per nucleon) at 
Tg = 9. In this study. Ye is taken as a free parameter 
and varied from 0.20 to 0.50 with an interval of 0.01. 

4.2. Nucleosynthetic Abundances 

Figure 4 shows the resulting neutron-to-seed ratio 
Yn/Yh at Tg — 2.5 (approximately the beginning of r- 
processing) for each (Mns, Li,) set as a function of Y^. 
Here, Y^ and 

Y,^ J2 ^(^'^) (9) 

Z>2,A 

are the abundances of free neutrons and the heavy nu- 
clei {Z > 2), respectively. The solid and dotted lines 
denote the results for Li/,51 = 10 and 1, respectively. For 
Mns = 1-4 Mq, Yn/Yh is only a few 10 at Y^ « 0.35 (that 
is the lowest Y valu e observed in hydrodynamic results, 
IWooslev et al.lll99l) . In this condition, oi ily the second 
r-process peak {A = 130) can be formed (jWanaio et al.l 

^ http:/ /nucastro.org/reaclib.html. We used the older version on 
the web, which has been recently updated with new experimental 
rates. 

^ http:/ /www. astro. ulb.ac.be/Html/bruslib. html. 



l2001f) . and a significantly lower Ye (~ 0.25) is needed for 
the third peak {A = 195) formation. On the other hand, 
for Mns = 2.0 M0, Y^/Yh ~ 100 is obtained with a rea- 
sonable choice of Ye (~ 0.4), whic h is sufficient for th e 
third peak {A = 195) formation (jWanajo et al.ll2001h . 
In the following, the results for Mns = 1-4 are pre- 
sented, keeping in mind that the values of Ye are sig- 
nificantly lowe r-shifted compared to the hydrod ynamic 
resuhs (Woosl ev et"allll994t lArcones et al.ll2007[ ). Note 
that qualitatively similar nucleosynthetic results can be 
obtained for Mns = 2.0 M0 with higher-shifted Ye val- 
ues. 

As seen in Figure 4, the Yn/Yh values in the early 
(^1^,51 — 10) and late (^1^,51 — 1) winds (in particular 
for Mns = 1.4 M0) at the same Ye are not significantly 
different, despite a large difference in entropy (Fig. 1). 
This is due to the shorter expansion timescale in the 
early wind, in which the entropy is lower. These two 
effects compensate with each other and lea d to similar 
nucleosynthetic results (IWanajo et al.l f2"001f ). We find, 
therefore, that the effect of Ye variation (that is most 
difficult to constrain) on nucleosynthesis is much larger 
than the Li, history. In fact, only ten percent change in 
Ye leads to a substantial difference in the r-process curve 
(Figs. 6 and 7). Note that the high Yn/Yh values at 
Ye = 0.48 - 0.50 for Mns = 2.0 Mq with L^,5i = 10 are 
due to its rather sho rt expansion timescale, in which Y/, 
is sig nificantly small (jHoffman et al.|[l997t IWanaio et al.l 

The time variations of r, Tg, p, and s for Rs — 215 km 
(= Rc), 500, 1000, 3000, 10,000, and 30,000 km in the 
early wind (-/ji/,51 = 10) with A/ns = 1.4 Mq are shown 
in Figure 5. The calculated abundance curves with these 
trajectories are displayed in Figures 6 and 7. The as- 
sumed Ye values of 0.35 ± 0.04 (Fig. 6) and 0.25 ± 0.04 
(Fig. 7) are responsible for the second and third peak 
formations, respectively. Such small variations in Ye 
with time would be inevitable when realistic time evo- 
lution of neutrino-d r iven outflows wer e considered (e.g., 
IWooslev et allll994t Uliras et al.ll2006f ). Roughly speak- 
ing, therefore, the envelope of these five curves in each 
panel might be regarded as the mass-averaged abundance 
pattern e xpected in mo re realistic, time-evolving winds 
see, e.g., IWanaiol [2OOI . 

In Figure 6, the envelope for Rs = 500 km best repro- 
duces the solar abundance pattern ( dots. [Kamelei et al.l 
Il989l ) between A — 100 and 180, including the second 
and rare-earth peak abundances. The Rs = Rc = 215 km 
case forms the troughs at the low {A — 110 — 120) and 
high {A = 140 — 150) sides of the second peak. This 
is in fact occasionally attributed to the p roblems in the 
predi cted nuclear masses around = 82 (jWanaio et al.l 
|2004| ). For Rs > 1000 km, the envelopes unacceptably 
broaden compared to the solar r-process curve. 

For the third peak abundances {A = 180 — 200 in 
Fig. 7), in contrast, the envelopes except for Rs = Rc = 
215 km are in reasonable agreement with the solar r- 
process curve. In particular, better fits are seen for 
Rs = 1000 km and Rs > 10, 000 km, while the intermedi- 
ate Rs (= 3000 km) case results in the lower-shifted third 
peak. Note that the abundance curves for Rs > 3000 km 
in Figure 6 and for Rg > 10, 000 km in Figure 7 cannot 
be distinguished. All these aspects can be understood as 
a consequence of the termination-shock deceleration of 



The r-Process and Termination Shock 



5 



supersonic winds, as discussed in § 4.3. 

4.3. Effect of Termination-Shock Deceleration 

The heavy nuclei synthesis ceases when the neutron- 
to-seed ratio decreases to Yn/Yh = 1 (hereafter, "n- 
exhaustion", IWanaio et all |2004 IWanaio .2007). Fig- 
ure 8 {top left) shows Tg at n-exhaustion, Tg^e, for each 
Rs as a function of Ye- In each line, a temperature jump 
appears as Ye decreases (except for Rs — Rc)- As an 
example, the line for Rs = 500 km {red line) overlaps 
with those for Rs > 500 km at Yg > 0.38. This is a 
consequence that n-exhaustion occurs before reaching Rs 
because of the low Yn/Yh values (Fig. 4). For Ye < 0.38, 
n-exhaustion takes place after the termination-shock de- 
celeration (with the temperature jump) owing to the suf- 
ficiently high Yn/Yh values. In these cases, Tg^c is almost 
independent of Ye because of the rather slow temperature 
decrease after the wind termination (Figs. 1 and 2). Note 
that, for Rg — 30, 000 km, n-exhaustion occurs before the 
termination-shock deceleration in all Ye cases. As seen in 
Figure 8 {top left), the Tg^c curves for Rg > 3000 km are 
overlapped at Ye > 0.30. Similarly, the Rs > 10,000 km 
cases result in the same Tg_e curve for Ye > 0.21. This 
is the reason why the abundance curves in these cases 
cannot be distinguished (Figs. 6 and 7). 

No heavier nuclei are synthesized after n-exhaustion, 
but local rearrangement of the abundance distribution 
continues until the neutron capture becomes slower than 
the competing /3-decay. The final r-proces s curve is thus 
fixed at Tn ^lrn = 1 (hereafter, "freezeoufi lWanajo et al.l 
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Here, 



^ Xn,{Z,A)Y{Z,A) 

J- h 



E 

Z>2,A 



and 



^ Xp{Z,A)Y{Z,A) 



Z>2,A 



(10) 



(11) 



are the abundance-averaged mean lifetimes of (n, 7) re- 
actions and /3-decays for Z > 2 nuclei, respectively, 
where Xnf{Z,A) and Xi3{Z,A) are the rates of the cor- 
responding reactions. These represent the lifetimes of 
the dominant species at a given time. The "classical" 
r-process is characterized by the conditions r„^/r-y„ w 1 
and t^/t„ ^ 1, where T^n is the (7,n) lifetimes defined 
similar to equation (11). 

The r-process path for each case can be specified in 
terms of the neutron separation energies. In this study, 
the abundance-averaged 2n separation energy S2n di- 
vided by two, defined by 



5a = 



1 



S2n{Z, A) 



Y{Z,A) 



(12) 



Yh ^ 2 

Z>2,A 

(|Wanaioll2007t l. is taken to represent the nucleosynthetic 
path at a given time. That is, S'a = is the neutron-drip 
line, while a higher S'a locates at closer to the /3-stability 
line. In Figure 8 {right; solid line), at freezeout, Sa.f, 
in each Rs case is shown as a function of Ye- The r- 
process paths predicted fro m the (n,7)-(7,n) approxi- 
mation (fCorielv fc Arnouldl[l996l ). 



S°(MeV) 



34.075 - \ogNn + ^ logTg j (13) 



are also shown in Figure 8 {right; dotted line) for com- 
parison purposes. 

As seen in Figure 8 {top right), the (n, 7)-(7, n) approx- 
imation can reasonably predict the r-process paths at the 
freezeout for Rs < 1000 km, in which 5° is systematically 
(0.5 — 1.0 MeV) higher than Sa,i- This is a consequence 
that the classical r-process conditions Tn-y/r-yn ~ 1 and 
Tp/Tn ^ 1 are approximately valid in these cases ow- 
ing to the sufficiently high temperatures and densities. 
In contrast, the (n,7)-(7,n) approximation poorly pre- 
dicts the real path at the freezeout for Rg > 3000 km. 
This is due to the low temperatures Tg_f < 0.7 as well 
as the low densities (Fig. 5) in these cases, in which the 
(n,7)-(7, n) e quilibrium is not valid ("cold r-process", 
IWanaiol I2007D . The classical conditions above are re- 
placed by Tn-y/Tjn «C 1 and T/3/r„ > 1. The latter con- 
dition, that is, the competing /3-decay with the neutron 
capture, pushes the r-process path S'a back to a higher 
value (see the similar situa tion in the decompos i tion of 
cold neutr on star matter, [L attimer ct al. 1977; Me veil 
17989 ; Freiburghaus et al.|[l999 ; Goriely ct al. 2005lJ). As 
a consequence, the Rg > 3000 km cases take similar 
Sa,f values. This is the reason that the Rg < 1000 km 
cases result in quite different abundance curves, while 
the other cases do not show significant differences. In 
addition, Sa.f values for Rg — 1000 km and > 10000 km 
are close at Ye ^ 0.25. As a result, the abundance 
curves around the third peak in these cases are similar 
(Fig. 7). Figure 8 {top panels) shows that the requisite 
conditions for the second and third peak abundances are 
Tg e = 1.2 — 1.4 and < 1.0, respectively, which result in 
Sa'f = 4.2 - 4.4 MeV and 2.8 - 3.8 MeV. 

In summary, in the early winds (Li/,51 — 10), the ter- 
mination shock can play a decisive role for the second 
and rare-earth peak abundances. Only a rather small Rs 
{^ 500 km) is allowed to reproduce the solar r-process 
curve in this atomic mass range {A « 100—180). The ter- 
mination shock may have a role for the third peak abun- 
dances as well (e.g., Rg = 1000 km). The solar r-process 
curve in this mass range {A « 180 — 200) is, however, well 
reproduced without shock (e.g., Rg = 30,000 km), too. 
Any subsonic winds cannot reproduce the solar r-process 
curve at all, in which n-exhaustion occurs at rather high 
temperature (Fig. 5; top right). 

In the late winds (L^ ^i — 1), all the nucleosynthetic 
results for Rg > Rc ~ 1056 km cannot be distingtushed 
(not displayed in this paper), which have similar abun- 
dance curves to those in the early winds (L^^si = 10) 
with Rg > 10,000 km (Figs. 6 and 7). This 'is due to 
the quite low temperature (Tg = 0.27, Fig. 1) at the 
sonic point. As a consequence, freezeout has already oc- 
curred at r < i?s in most cases (Fig. 8; bottom left), 
in which the termination-shock has no effect on the r- 
process. This leads to almost the same Sa,f curves (Fig. 8; 
bottom right), which are similar to those in the early wind 
(L^,5i ^ 10) with Rg > 10,000 km (Fig. 8; top right). 
The resuhing Saj values (2.9 - 3.9 MeV) faU within the 
requisite range to reproduce the solar r-process curve 
near the third peak {A « 180 — 200), but too small for the 
second and rare-earth peaks. Therefore, the late winds 
can be responsible only for the third peak abundances 
and heavier, in which the effect of termination shock is 
of no importance. It should be noted that subsonic winds 
may be able to reproduce the second and rare-earth peak 
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abundances, if Tg^c appears to be 1.2 — 1.4 (Fig. 1; top 
right). 

4.4. Effect of Termination- Shock Heating 

As noted in § 3, the large entropy jumps owing to 
the termination-shock heating do not help to increase 
Yn/Yh, since the jumps appear only when the tempera- 
ture decreases below Tg — 2.5 (Fig. 3). In fact, all the Rs 
cases for a given (A/ns, L^) set have the identical Yn/Yh 
curve as a function of Ye in Figure 4. One may con- 
sider, however, a large entropy jump during r-processing 
affects the nucleosynthetic path and modifies the final 
abundance curve. An entropy increase is equivalent to 
a reduction of the matter density (and thus the neutron 
number density) at a fixed temperature in the current 
radiation dominated outflows. A sizable entropy jump 
can thus, in principle, modify the r-process path. 

Surprisingly, even the entropy jump to about 400 N^k 
does not change the final abundance curves. In Figures 6 
and 7, the dotted lines denote the nucleosynthetic results 
by suppressing the entropy jumps without changing the 
temperature histories (but artificially increase the densi- 
ties). No differences between the abundance curves with 
and without entropy jumps can be seen in Figure 6, and 
only minor differences appear for Rg = 3000 km and 
10,000 km in Figure 7. 

The reason is that the entropy jump is too small in 
the Rs < 1000 km cases (Fig. 3) to change the r-process 
path during the freezeout phase. On the other hand, a 
large entropy jump appears only at late times when the 
r-process has already ceased in most cases. As an exam- 
ple, for Rs — 30,000 km, all the Ye > 0.20 cases experi- 
ence n-exhaustion before the termination-shock heating 
(Fig. 8, top left). For Rs = 3000 km and 10,000 km, the 
winds with < 0.29 and Y^ < 0.21, respectively, reach 
the termination shock radii before n-exhaustion. These 
cases, however, result in only minor modifications of the 
final abundance curves (Fig. 7). 

5. CONCLUSIONS AND DISCUSSION 

We have investigated the role of wind termination 
shock on the r-process nucleosynthesis in supersonically 
expanding neutrino-driven outflows. For this purpose, 
a semi-analytic, general relativistic neutrino-driven wind 
model was utilized to obtain spherically symmetric tran- 
sonic wind solutions. The Rankine-Hugoniot relations 
were applied at arbitrary chosen shock radii Rs to obtain 
the boundary conditions for the subsequent termination- 
shocked subsonic outflows. 

We explored the properties of termination-shocked su- 
personic winds for various neutron star masses Mns, neu- 
trino luminosities L^, and termination-shock radii Rs- 
Mns was taken to be 1.4 M© and 2.0 Mq as a typical 
and a most massive proto- neutron stars, respectively. 
Lj^^5i was taken to be 10 and 1 as representative of early 
and late winds, respectively. The earlier winds with 
Lu.5i > lOi iti whi c h matter is likely to be proton-rich 
(iBuras et al.l l2006t iFrohlich et all l200d iKitaura et aH 
|2006[) . mav not be relevant for the r-process. The later 
winds with L^^^i < 1 may not contribute to the Galac- 
tic production of r-process nuclei owing to their small 
mass ejection rates ( < lO'^M©. IWanaio et al.l [20011 : 
[Thompson et al.ll2001[ ). 



5.1. Conclusions 

Rs was changed from the sonic radius i?c 200 — 
1000 km) up to 30,000 km. Our main results are sum- 
marized as follows. 

1. The shock jumps are greater at larger Rs from the 
sonic point, for fixed Mns and L^. The termination- 
shocked winds become, by definition, subsonic outfiows, 
in which the temperature and density decrease rather 
slowly with radius (and with time owing to the decay- 
ing expansion velocity). The entropy increases by the 
termination-shock heating up to several 100 NAk in some 
extreme cases (Rs = 30,000 km), which is several times 
larger than its unshocked value. The jump is, on the 
other hand, only moderate at small Rs (^ 1000 km), 
which disappears in the (subsonic) Rs = Rc case. 

2. For fixed Mns and Rs, the shock jump is more 
prominent in the earlier wind (ii/,51 = 10) owing to 
its faster wind velocity. In the early wind, the shock 
jump occurrs even at Rs — several 100 km when the 
temperature is as high as Tg > 1, because of its smaller 
sonic radius {Rc ~ 200 km). On the other hand, the 
jump occurrs only at low temperature (Tg < 0.2) in the 
late wind {L^^^i = 1) owing to its distant sonic radius 
(Rc ~ 1000 km). For fixed and Rs (as well as for a 
fixed i?NS = 10 km), the jump is significantly greater in 
the massive Mns (= 2.0 Mq) case owing to its somewhat 
faster wind velocity. 

3. The entropy jump takes place only after the tem- 
perature decreases below Tg « 2.5, except for extremely 
high L^. Therefore, the termination-shock has no effect 
on increasing the neutron-to-seed ratio. 

All these results are consistent, at least qualitatively, 
with those found in one-dimensional hydrodynamic sim- 
ulations by Arcones et al. (2007). We further performed 
detailed nucleosynthesis calculations to examine in par- 
ticular the effects of shock deceleration and shock heating 
on the r-process. 

4. The temperature in the shock decelerated wind de- 
creases rather slowly with time. As a result, the tempera- 
ture stays mostly con stant during the r-process freezeou t 
phase (as assumed in lWanaio et al.|[2002l : IWanaioll2007D . 
We find that the decelerated temperature is highly de- 
pendent on Li, and Rs. In the early wind (Ti/,51 = 10), 
the termination-shock deceleration plays a decisive role 
to determine the final r-process curve. A solar-like r- 
process curve including the second and rare-earth peaks 
(A « 100 — 180) can be obtained with rather small Rs 
(~ 500 km), in which the r-process freezeout takes place 
at Tg « 1.2 — 1.4. On the other hand, the heaviest por- 
tion of the solar r-process curve {A « 180 — 200) can 
be reproduced both with Rs ^ 1000 km (Tg w 1.0) and 
Rs > 10,000 km (Tg < 0.3). In the latter cases, the 
r-process has ceased before the shock deceleration, in 
which the termination shock plays no role. In the late 
winds (Ti,^5i — 1), the shock deceleration can occur only 
at Rs > 1000 km, in which the r-process has already 
frozen out in most cases. In the late winds, therefore, 
only the heaviest part of the solar r-process curve can 
be reproduced, in which the termination shock has little 
effect. 

5. The entropy jumps by the termination-shock heat- 
ing have negligible effects on the final r-process curve. 
This is a consequence that a sizable entropy jump occurs 
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only for Rg > 10000 km when the temperature is as low 
as Tg < 0.3, where the r-process has already ceased. For 
smaller Rg, the shock heating can occur before freezeout, 
but the entropy jump is too small to modify the r-process 
path. 

5.2. Discussion 

Our approach was based on the simple steady state 
wind solutions with arbitrary chosen (but reasonable) 
input parameters, although our results are qualitatively 
consistent with mor e detailed hydrodynamic results by 
lArcones et al.l (|2007D . We also employed rather low ini- 
tial Ye (= 0.2 — 0.4) to obtain the r-process abundances 
with reasonable parameter settings. Our results pre- 
sented here should be therefore taken to be only sugges- 
tive. Keeping this caution in mind, we suggest a couple 
of possible supernova progenitors that may be relevant 
for the astrophysical r-process site. 

First is a rather massive progenitor 25 Mq), in 
which the termination-shock radius may s t ay cfo se to 
the sonic point for long time. lArcones et al.l ()2007f ) show 
that a more massive progenitor has a higher mass ac- 
cretion rate that results in a higher and a smaller 
Rs- As a consequence, in their 25 Mq case, the tem- 
perature stayed at Tg > 1.0 during the first 10 seconds 
after core bounce. In this condition, both the lighter 
(A « 100 - 180) and heavier {A « 180 - 200) parts 
of the solar r-process curve would be reproduced if req- 
uisite neutron-to-seed ratios were achieved. This con- 
dition Tg > 1.0 (as well as high neutron number den- 
sity, > 10^° cm~'^) has in fact long been considered 
as a physical requirement to account for the solar r- 
process abundanc e curve fe.g.. lMathews fc Cowa^ll990l : 
iKratz et al.lll993l ). A massive progenitor forms a massive 
proto-neutron star, in which general relativistic effects 
may in part help to in crease ent ropy an d reduce expan- 
sion timescale fcardall fc Fulleiiri 997: Ot suki et al.ll2000l : 
IWanaio et aI1l2001l : IThompson et al.ii200]E 

Second is, in contrast, a star of ^ 10 M p, near the low- 
mass end of supernova-progenitor range. lArcones et al.l 
()2007t ) show that, in their 10.2 M© case, the termina- 
tion shock propagated outward quickly owing to the 
steep density gradient, in which the temperature at Rg 
dropped to Tg ^ 0.1 after 10 seconds. In this case, 
the lighter part {A « 100 — 180) of the solar r-process 
curve should be formed in shock-decelerated early winds, 
in which the (shocked) temperature is still as high as 
Tg ss 1.2 — 1.4. The heavier part may be followed 
in later winds, in which the r-pro cess proceeds in very 
low temperatures (cold r-process, iWanaiol [20071) . Low- 



mass supernovae have been also suggested as the as- 
trophysical r-process site, fro m Galactic chemical evo- 
lution studies (~ 8 - 10 Mr;,, iMathews fc CowanI 119901: 
llshimaru fc Wanaiol Il999t llshimaru et all l2004f ). The 
mechanism to obtain requisite physical conditions for the 
r-process in such a supernova, however, has been lacking. 

Needless to say, we cannot exactly constrain the mass 
range of supernova progenitors on the basis of our cur- 
rent (only suggestive) results. It is likely, however, that 
the (shocked) temperatures in intermediate progenitor 
cases (e.g., 15 - 20 M©) stay at Tg ~ 0.5 - 0.7 for long 
time, in which the third peak abundances may be lower- 
shifted compared to the solar r-process curve (Fig. 7, 
see also IWanaio et al.ll2002tlWanaioll2007l ). It should be 
also noted that the subsonic outflow has higher temper- 
ature than that in the transonic wind at any radii. We 
cannot exclude, therefore, possibilities that these progen- 
itors provide s uitable temperature Tg « 1.0 — 1.4 during 
the r-process. lArcones et al.l (|2007( l show, however, that 
all their simulations with "standard" parameter settings 
resulted in forming supersonic outflows within 10 sec- 
onds after core bounce. Only a subset of simulations 
with somewhat "extreme" parameter choices appeared 
to have the transiting features between transonic and 
subsonic outflows. 

Given that our speculations above are correct, it is 
still difficult to specify which of ~ 10 Af© or ~ 25 M© 
stars are the dominant sources of the Galactic r-process 
nuclei. Both cases can potentially reproduce the whole 
mass range of the solar r-process curve (but Pb abun- 
dance can be substantially different, see Fig. 7 and 
IWanaioll2007f ). Detailed nucleosynthesis calculations in 
time-evolving neutrino-driven winds with the termina- 
tion shock will be needed to draw more definitive an- 
swers. It should be also noted that we assumed only 
spherically symmetric outflows in this study, while re- 
cent hydrodynamic studies suggest some decisive roles of 
multidimensional effect s on the mechanism of neutrino- 
driven explosions (e.g.jBuras et alll2006t lOhnishi et al.l 
l2006t iBurrows et al.l 2007h . Eventually, we will need ex- 
tensive nucleosynthesis studies in the outflows obtained 
from multidimensional core-collapse simulations. 
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Fig. 1. — Radial velocity u {top left), temperature T {top right), density p {bottom left), and entropy s {bottom right) in neutrino-driven 
winds for Mns = 1-4 Mq as functions of radius. The neutrino luminosities are taken to be L^^5i = 10 {blue) and 1 {green), respectively. 
In each Li, case, the termination-shock radius is varied from the sonic radius (denoted by star) to infinity (see text). 
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Fig. 2. — Same as Figure 1, but for Mns = 1-4 Mq (green) and 2.0 Mq (red). The neutrino luminosity is taken to be L^,5i = 10. 
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Fig. 3. — Relations between the entropies Ss and the temperatures Tg just above the termination-shock radii Rb- The neutron star masses 
and the neutrino luminosities are taken to be M^s = 1.4 Mq (crosses), 2.0 Mq (asterisks) and i/i/,6i = 10 (solid lines), 1 (dash-dotted 
lines), respectively. For A/ns = IAMq, Lv,5i = 70 (dashed line) is also considered as an extreme case. Rs in each Lv case is varied from 
Rc (sonic radius, marked by circles) to 500 km, 1000 km (marked by squares), 3000 km, 10000 km (marked by triangles), and 30000 km, 
whenever Rs > Rc is satisfied. 




Fig. 4. — Neutron-to-seed ratios Yn/Yf^ at the beginning of r-processing (Tg = 2.5) in Mns = 1.4 Mq and 2.0 Mq cases as functions of 
the initial electron fraction Yg. Solid and dotted Unes denote the results for 51 = 10 and 1, respectively. 




Fig. 5. — Time variations (t = at Tg = 9) of r {top left), Tg {top right), p {bottom left), and s {bottom right) for the selected shock 
radii Rs- The neutron star mass and the neutrino luminosity are taken to be Mns = 1-4 Mq and = 10, respectively. 
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Fig. 6.— Comparison of the nucleosynthetic yields (lin es) for Ys = 0.31, 0.33, 0.35, 0.37, and 0.39 in the (Mns.^i^.Si) = (1.4M0,1O) 
case with the solar r-process abundances f dofs: lKappeler et al.lll989l) as functions of mass number. The latter is shifted to match the height 
of the second peak (A » 130). 
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Fig. 8. — Temperatures at n-exhaustion Tg^c {left) and nucleosynthetic paths at freezeout (in terms of neutron separation energies) Sa,f 
(right) for Mns = IAMq as functions of Ye. The termination-shock radius Rb is varied from the sonic radius up to 30000 km. Top and 
bottom panels show the results for L^,5i = 10 and 1, respectively. The paths predicted from the (n,7)-(7, n) equilibrium are also shown 
by dotted lines (right panels). 



